Robert Warneford-Thomson
pacman::p_load(
colorspace,
ggbeeswarm,
ggnewscale,
ggh4x,
ggpubr,
ggsci,
magrittr,
paletteer,
patchwork,
scico,
tidyverse
)
knitr::opts_chunk$set(
echo = FALSE,
fig.path = "./fig-1/",
fig.dpi=600,
warning = FALSE,
message = FALSE,
dev = c("png", "pdf")
)
setwd("~/Analyses/bms-dms/paper")
cbPalette <- c("#DC5E65", "darkgrey", "#0072B2", "#56B4E9", "#E69F00", "#009E73", "#F0E442", "pink", "#CC79A7", "lightgrey", "grey")
aa_order <- c(
"*", "P", "G", "A", "M", "V", "L", "I", "T", "S",
"C", "Q", "N", "Y", "W", "F", "E", "D", "K", "H", "R"
)
source("../../dms/src/model/dms-analysis-utils.R")
After discussing with Conor, will use 1 and 10 U/mL IFNa data from run3, and 100 U/mL IFNa from run7
Iām going to just use a 1% FDR to determine GoF and LoF variants with no minimum log2fold or midpoint shift threshold
FigureĀ 1